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Abstract 



For time projection chambers the accuracy in measurement of track coordinates along the pad-row 
direction deteriorates with the drift distance (z): cxx ~ D^ ■ z/Nes, where D is the diffusion constant 
and A'^eff is the effective number of electrons. Experimentally it has been shown that A^eff is smaller 
1-^ ' than the average number of drift electrons per pad row {{N)). In the previous work we estimated 
J^ . iVeff by means of a simple numerical simulation, taking into account the diffusion of electrons only 
Ph ! in the pad-row direction [l[. The simulation showed that iVefj could be as small as ~ 30% of (A^) 
^ ' because of the combined effect of statistical fluctuations in the number of drift electrons (A^) and in 
. ^ ' their multiplication in avalanches. In this paper, we evaluate the influence of the diffusion along the 
£^ ' direction of pad-row normal on the effective number of electrons. 
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1. Introduction 



In the previous paper we estimated the effective number of electrons (A^eff) contributing to the 
coordinate measurement of a time projection chamber (TPC) equipped with a Micro-Pattern GAS 
cr^ • Detector (MPGD) and ideal readout electronics [l]. A"eff parametrizes the spatial resolution for a 
^~~' . pad row as follows 

where ax is the spatial resolution along the pad-row direction, axo is the intrinsic resolution, and D 
denotes the transverse diffusion constant, with z being the drift distance. In the ideal case (with an 
infinitesimal pad pitch) o"xo vanishes for particle tracks perpendicular to the pad row, and increases 
with the track angle (0) with respect to the pad-row normal because of the angular pad effect. Only 
right angle tracks (0 = 0°) are considered throughout the present work. 
Under the conditions listed in Ref. [l[, A'eff is given by 

^ /iV(l + /) (2) 



A^eff \N/ 

Iril 
with q being the signal charge induced on the pad row by single drift electron^j. Although Eq. 



where A^ denotes the total number of drift electrons detected by the pad row and / = cr"^/ {q) , 



1 It should be noted that Eq. © is an expression for 1/N{h) in Eq. (7.33) of Ref. [i|. 
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was derived assuming the total charge (J2f=iQi) to be constant (A^ ■ (g)) it was found to be a good 
approximation by a numerical simulation for a practical value of (A^) (see Ref . pj] , and Appendix A 
for details). 

The simulation for argon-based gas mixtures gave A^cff of ~ 22 for / = 2/3 and a pad-row pitch 
of 6.3 mm [l|], which is about 30% of the average value of A^ {{N) ~ 71), and is consistent with 
the values obtained with a small prototype TPC jsMsj . The value of Nes corresponds to ~ 36 for a 
pad-row pitch normalized to 1 cm, assuming N^g to be (approximately) proportional to the pad-row 
pitchl. 

Recent resolution measurements with a larger prototype TPC with MicroMEGAS readout, how- 
ever, gave a significantly larger estimate for N^s (~ 56 for 1-cm pad height) [6jj. A possible origin 
of the discrepancy could be the de-clustering of drift electrons due to diffusion along the direction 
of the pad- row normal (-Dy), which is expected to be more efficient for larger TPCs with a longer 
average drift distance. It should be pointed out that the diffusion only along the pad-row direction 
(Dx) was taken into account in Ref. [l|. In the present work, we evaluate the contribution of the 
de-clustering effect to the increase of A^eff, through the decrease of {l/N) due to the finite Dy. 

A hypothetical case of an infinitely large TPC is considered in Section 2, a more practical case is 
qualitatively discussed in Section 3, the results of a realistic simulation are shown in Section 4, and 
Section 5 concludes the paper. Readers are suggested to have read Ref. [l| in advance. 

2. In the case of infinitely large TPC 

Let us consider the hypothetical case of a TPC with infinitely large drift distance and readout 
plane. Primary ionization clusters created along a particle track at infinitely large drift distance get 
completely de-clustered and the total number of drift electrons reaching a pad row obeys Poisson 
statistics with a mean fi = (N): 

The dimensions of the readout pad plane are considered to be large enough. Otherwise a part of the 
drift electrons created at long drift distances would be absorbed by the field cage (the inner or outer 
wall of cylindrical TPC) before reaching the readout plane. 

The average value of the inverse of A^ in this case is given by 

1 \ °° 1 

n) = S^-W 

^^ ' N=l ^^ 



N ■ N^ 



= (E,(/x)-ln(/i)-7)-e-^ 

where Ei is the exponential integrajjand 7 is the Euler-Mascheroni constant (~ 0.577). When A^ = 
the pad row is inefficient and provides no coordinate measurement. Therefore it is excluded from 



^ Actually, this is a bold assumption. See Appendix A. 

•^ In fact, the values of N^s were obtained using different kinds of charged particle: a beam of 5-GeV/c electrons 
in Ref. [6] while a beam of 4-GeV/c pions or cosmic rays in Refs. [3|-|5|- The discrepancy is still large, however, even 
if the difference in the primary ionization density is taken into account. 

"* The exponential integral is defined as 

r e* 

E^{x) = -dt. 

J —00 ^ 



the summation. Fig. [T] shows the behavior of i?N = (N) ■ (l/N) as a function of (A^). For a practical 
value of (A^) > 50, {1/N) ~ 1/ (A^) and A^eff is expected to be ~ (A^) /(I + /). 




Figure 1: R^ = (TV) • (1/-/V) as a function of the average number of electrons in the case of Poisson distribution for N. 

For a finite drift distance A^eff < (N) /(I + /) because of broad and asymmetric distribution of 
A^ (Landau tail) even for relatively large (A^). For long drift distances the distribution approaches 
a Poissonian, which is close to a Gaussian for (A^) > 20. Consequently, A'efr is expected to be an 
increasing function of the drift distance [z) with an asymptotic maximum of ~ (A^) /(I + /). 



3. In the case of real TPC 

We evaluate in this section the variance (Var) of the total number of electrons reaching a readout 
pad row (A^) since Var(A^) gives a good measure for the deviation of A^-distribution from a Poissonian, 
for which Var(A^) = A^. First, let us suppose that a single point-like electron cluster of size n is 
created at a coordinate y = Y in the direction of the pad-row normal and z = Z in the drift direction, 
measured from the readout plane. The electrons diffuse on their way towards the readout plane. Their 
spread in the y-coordinate is given by a Gaussian with a standard deviation oi a = ay = D ■ y/Z, at 
the pad rows with a height of /tl- The probability to find u electrons reaching the pad row is given 
bjS 

in 



P{u-Y,a) 



n'^(r) ■ (1 - n(F))' 



where 



with 



h/2 

U{Y)=I G{y;Y,a)dy 

-h/2 



G{y;Y,a) = 



27Ta 



■ exp 



2a2 



^ More precisely h should be understood as the pad-row pitch, which is usually slightly larger than the pad height 
when the readout plane is covered over with pads. The pad-row pitch and the pad height {h) are not distinguished in 
the present paper. 

^ Hereafter the notation /"(x) represents the n-th power of a function /(x), i.e. (/(a;))". 



Since P(i^; Y, a) represents a binomial statistics for a fixed Y (and a), {y) and (z/^) are given by 

(z/(F)) =n-n(F) 

and 

^z/2(y)^ = n ■ n(F) + n(n - 1) ■ T{^{Y) . 

Let us further assume that the initial cluster is randomly created in a |/- region [—H/2, +H/2] 
with if > cr (see Fig. [2]). Then, averaging over -ff/2 <Y< +H/2, 



y = Y 



Pad row is assumed to be 
aligned with the X axis. 



Electron distribution 
after diffusion 




-H/2 



-h/2 +h/2 +H/2 



Figure 2: Illustration of the relevant variables. 
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where 



^\hj \2a 0F h 



( h 
exp 



4(t2 



See Appendix B for the derivation of function g{a/h). 

If there are A^cl (independent) clusters in [—H/2, +H/2], the average and the variance are given 
by multiplying A^cllJ: 

Taking the limit of if — ?■ oo while keeping the cluster density p = Ncl/H constant, 

(u) = p ■ n ■ h 
Var(z/) = p ■ n ■ h + p ■ n{n — 1) ■ h ■ g { — \ . 

In reality the cluster density (p) depends on the cluster size (n): 

p = p{n) = po-p{n) 

where po is the total primary ionization density and p{n) is the proportion of the cluster of size n 
{J2nP{''^) = !)• The number of electrons detected by the pad row (A^) is the sum of the contribution 
(i/) from various cluster sizes. Consequently, its average and variance are given by 

(A^) = po- h-'^n- p{n) 



Var(A^) = pQ-h-Y,n-p{n)-[l + {n~l)-g(^ 
Fig. [3] (a) shows 



Var(z/) ^ , ^, fa 

^ - 1 + (^ - 1) ■ ^ r 



for several values of the cluster size n as a function of the scaling parameter a /h. It is clear that 
V distribution approaches a Poissonian (Var(z/)/ (z/) = 1) slowly with the increase of a /h, especially 
for large clusters. Fig. [3] (b) shows the variance divided by the average for a realistic A^-distribution: 

Var(Ar) _ Unn ■ p{n) ■ {l + {n - I) ■ g (ff) 
(N) J2n n ■ pin) 

^ l^ En^(^-l) -pin) _ /a 
T,nn-p{n) \h 

along with the ratios obtained with a numerical simulation (see below). The probability mass function 
p(n) was assumed to be that corresponding to the cluster-size distribution shown in Fig. 2 of Ref. p|. 



Note that 



/WCL \ 

^ Mi \ = NcL ■ {v) and 
\ j=i / 

/JVCL /JVCL \ \ ^\ / /Nc 



£^'-(£^0 ) = ( E('^'-(^»)) ) = ^cL-((^-M)') 



vi— 1 \i— 1 // / \\i^l 
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Figure 3: Variance divided by mean as a function of the scaling parameter a/h: (a) for several fixed values of the 
cluster size and (b) for a realistic cluster-size distribution along with the data points given by a numerical simulation 
(filled circles). 



4. Evaluation of Nefi by simulation 

The analytic approach described in Sections shows that N^s is expected to be a slowly increasing 
function of cr/h, i.e. the drift distance. In order to confirm this quantitatively, we evaluated A^eff by 
means of a numerical simulation. 

The simulation code is identical to that used in the previous work pL|], except that the diffusion 
of drift electrons in the the direction of pad-row normal (-Dy) is taken into account. Initial electron 
clusters are randomly generated along the y-axis (with the pad row aligned with the x-axis) in a 
range wide enough compared to the diffusion {a = D ■ vZ) and the pad height (h). The cluster 
density is assumed to be 24.3 cm~^ x 1.2 (relativistic rise factor) as in the previous paper [l|. The 
size of each cluster is determined randomly using the probability mass function p{n) (see Section 
3). The secondary electrons originated from each cluster are then dispersed in the directions of 
pad row (x) and pad-row normal (y) with ax = Cy = a = D ■ yZ |3. The electrons with the final 
position located within the pad row (\y\ < h/2) are accepted (see Fig [2]). Gas gain is assigned to 
each accepted electron randomly assuming a Polya distribution {6 = 0.5, corresponding to / = 2/3) 
for the avalanche fluctuation. The coordinate resolution (o"x) is calculated using the charge centroid 
of the accepted electrons in the pad-row direction (x). The square of the ratio of the diffusion (a) 
to the resolution (o"x) gives Nes from Eq. ([1]) with crxo = 0. 

Fig.Hlshows the obtained Nqs and {1/N)~ as a function of a/h. The effective number of electrons 
certainly increases with a/h in association with the increase of {1/N)~ . However, the increase of 
iVeff is rather slow and would be observable only for large values of a/h, i.e. at (very) long drift 
distances. Examples of the resolution squared as a function of the drift distance are shown in Fig. \5\ 
for pad heights of 6.3 mm and 3.15 mm. The chamber gas is taken to be Ar-CF4 (3%)-isobutane 
(2%) as in the experiments of Refs. |5|, Isj]. The deviation from the linear dependence ( Eq. ([1]) with a 
constant N^s) is prominent at long drift distances without magnetic field, in particular for a shorter 
pad height. 



Note that D^ ~ Dy = D since the magnetic field (if it exists) is parallel to the z-axis. 




Figure 4: Inverse of the average of 1/A^ and effective number of electrons as a function of a/h for h = 6.30 mm and 
3.15 mm. The values for small (realistic) a/h are shown in (b). 
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Figure 5: Resolution squared as a function of the drift distance (z) given by the simulation for B = and 1 T, with 
a pad height (h) of 6.3 mm. The diffusion constant is assumed to be 315 (101) /ini/-y/cm for B= (1) T given by 
Magboltz [7[ for a gas mixture of Ar-CF4 (3%)-isobutane (2%) and a drift field of 250 V/cm. The relation simulated 
for B = T with h = 3.15 mm is also shown for comparison. The dashed straight lines show the linear relations with 
the values of NoB fixed to those at z = 0, in the case of B = T. For i? = 1 T, iVoff is almost constant throughout 
the drift distance shown in the figure. 



5. Conclusion 

We estimated the effect of the electron diffusion along the direction of pad- row normal (-Dy) on 
the spatial resolution of TPCs operated in argon-based gases. It does affect the effective number of 
electrons contributing to the azimuthal coordinate measurement. The value of N^s increases with 
the drift distance since the distribution of the number of electrons detected by a pad row approaches 
a Poissonian by de-clustering. However, the de-clustering process is rather slow because its scaling 
parameter is cr/h, and N^ff can be assumed to be constant for practical TPCs. 

Furthermore, the effect of avalanche fluctuation (Rq) was found to be almost constant (~ 1 + /) 
for a realistic pad-row pitch greater than ~ 6 mm (see Appendix A). Therefore Eq. ([2]) is expected 
to be an appropriate expression for the value of iVefj, with (l/N) estimated with Dy = 0. 

It is unlikely that the increase of A'^efr observed with a large TPC prototype, with the maximum 
drift length of ~ 600 mm, arises from the finite Dy. The larger N^s may have been owing to other 
factors such as smaller avalanche fluctuations (/), or improvement of the signal-to- noise ratio and/or 
better calibration of the readout electronics (see Appendix C of Ref. |5|). 

The increase of Nes would be observed at long drift distances with a large TPC operated in a 
gas with a relatively large transverse diffusion constant in the absence of magnetic field. 
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Appendix A. Pad-height dependence of Nef[ 

In the previous work, the pad-row pitch (~ pad height h) was fixed to 6.3 mm [1]. We show here 
the pad-height dependence of the effective number of electrons given by a numerical simulation. The 
simulation code is exactly the same as that developed for Ref. |ll]. Therefore the diffusion of drift 
electrons in the direction of pad-row normal (-Dy) is not taken into account and the values of Ng^ 
are those for zero drift distance. 

If we write A^eff = (A^) /R, with R being a reduction factor, R is expressed as R = -Rn ■ Rq, where 
i?N = (A^) ■ {N~^) and -Rq derives from the avalanche fluctuation in the detection device for each 
of drift electrons. The value of -Rq is expected to be close to 1 -|- / for large (A^) (large pad height) 
since I]^i g, ~ A^ ■ (q) becomes a good approximation. Fig. lA.ll shows the pad-height dependences 
of -Rn, -Rq and -R. The relative variance of avalanche fluctuation (/) is assumed to be 2/3. The value 
of -Rq is almost constant (~ 1 + /) for practical pad heights ( ^ 6 mm) whereas -Rn is a decreasing 
function of the pad height as expected. 

The values of {1/N)~ (= (A^) /-Rn) and Aefr are plotted in Fig. lA. 21 against the pad height, along 
with (A^). It is clear that Agfr is not a linear function of the pad height because of -Rn decreasing 
with the pad height. The effective number of electrons is about 31% (43%) of (.A) for a pad height 
of 6.3 mm (100.8 mm). 
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Figure A.l: Reduction factors {R, Rn and Rq) as a function of the pad height. 
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Figure A. 2: (a) Average number of electrons {{N)), inverse of the average of 1/A^ {{1/N) ), and effective number of 
electrons (A'cff ) as a function of the pad height, (b) Details for small pad heights. 



Appendix B. Derivation of function g 



In this appendix we derive the explicit expression of the function g{a/h) used to evaluate Var(z/) 
in Section 3. The function is defined as 



/(j\ rH/2 roo 

h-gl-)= lira U^(Y) dY = U^(Y) dV 

\hj H^ooJ-H/2 J-oo 



where 



with 



U{Y) 
G{y;Y,a)^ 



h/2 
-h/2 



G{y; Y, a) dy 



2na 



■ exp 



{y - yf 

2(t2 



Let us carry out the integration on the right hand side of the equation: 



n2(r) dY 



Itxcj'^ 



2^a'^ 



oo 
oo 



h/2 
-h/2 



-4-^-^h-i:M-^^h' 



dY 



/■oo rh/2 rh/2 f 1 / \ 

dY dy dyf exp -— ({y - Yf + {yl - Yf) 

J-oo J-h/2 J-h/2 . ZO-^ ^ ^ 



-h/2 J-h/2 
rh/2 ph/2 



I rOD j'h/2 rh/2 

— - / dY \ dy dyf exp 

ra'^ J-oo J-h/2 J-h/2 



27ra- J-oo 
1 /-Vs 



'h/2 J-h/2 
rh/2 



2^2 
1 



y + yfy 1 2' 



dy \ dyl exp 

-hl2 J-h/2 



2y^Cr J-h/2 ' J-h/2 
2 /-/i/v^ fh/y/2-i. 



r I ^^ 

V7ra JO Jo 



4^2 

2 N 



exp 
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2a2 



{y - yf 



dr] 



with ^ = -= — and r] — 



2h /•! , fh-{l-u)/2a 

—= du exp(— w ) dv 

'IT Jo Jo 



V2 



with u = ——C and v = _ 
h ^ y/2a 



V2 



h / f{u) du 
Jo 



with f{u) 



2 rh-{l-u)/2a 



TT Jo 



exp(— f ^) dv 



h- ([m/(m)]J- / uff{u) du 



h^ /-i 
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u ■ exp 
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2^ 
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l.h/2a / 2(7 \ , 2N , 
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/ t ■ exp(-t2) rft 



10 



= h- 


eri — ^- / expf 

\^ \2a ) y/irh Jo 
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—w) dw 
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= t' 


= h- 


\ J h\ 2a f ( h^\ M 
\2a J Vtt/i \ \ 4cr^/ / 




Hence, 




fa\ r(h\ 1 2a 
\hj \2o"y V^ h 


1 — cxp 




It should be noted that the error function is given by 




erf(x)= ^^ re-'' dt 

VTT JO 




and that 









(?(0) = erf (oo) = 1 and g{po) = . 
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